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Transport properties, spectral function and optical conductivity of the adiabatic one-dimensional 
Su-Schrieffer-Heeger (SSH) model are studied with particular emphasis on the model parameters 
suitable for Rubrene single crystals based field effect transistors. We show that the mobility, cal- 
culated by using the Kubo formula for conductivity, vanishes unless we introduce an "ad hoc" 
broadening of the system energy levels. Furthermore, the apparent contradiction between angle 
resolved photoemission data and transport properties is clarified by studying the behavior of the 
spectral function. Finally, a peak in the optical conductivity at very low energy is obtained and 
discussed in connection with the available experimental data for Rubrene based devices. 



I. INTRODUCTION 

The organic ficld-cffcct transistors (OFET) play an im- 
portant role in the field of the so-called plastic electronics. 
Recently, OFET, based on single crystals of oligoacene 
molecules, have been developed^ From the study of these 
systems it stems that charge mobility in oligoacene sin- 
gle crystals can be one order of magnitude larger than 
that one with non ordered molecules. Among them the 
more promising are those based on Rubrene crystals that 
exhibit a strong anysotropy and the largest mobility mea- 
sured in organic semiconductors. 

In spite of the important technological impact of such 
devices, the intrinsic transport mechanism acting in 
rubrene is not fully understood. Actually, while the tem- 
perature dependence of the mobility exhibits a power- 
law (fi e ~ T~ s with 8 ~ 2)2r— reminding that one of 
mobile charge carriers (band transport), computational 
data indicate that the scattering length of the charge 
carrier becomes too short to be compatible with band 
transport. 5 Furthermore, some spectroscopic evidences 
support the localization of the charge carrier within one 
or few molecules^ and, on the other side, angle resolved 
photoemission spectroscopy (ARPES)£ seems to point 
again towards a band-like behavior of the charge carriers. 

In this framework a very interesting and simple one- 
dimensional model has been recently introduced by 
Troisi 8 where the charge carriers interact with the inter- 
molecular modes leading to a modulation of the charge 
carrier hopping. The proposed model is somehow very 
close to the SSH model introduced in a different context. 9 

The model has been studied within the adiabatic limit 
(phonons obey a classical dynamics) in Ref.Q by us- 
ing an approximated dynamical approach and the results 
have been, then, confirmed and extended by Ciuchi and 
FratiniiS who used, instead, a thermodynamic approach 
where the vertex renormalizations are neglected. The 
main results claimed by those authors is that a power 
law for the mobility temperature dependence can be re- 
covered within the proposed model and that the charge 
carriers involved in the transport undergo a "dynamical 
localization" . 

From this brief discussion it is clear that a systematic 



study, "numerically" exact, of the transport properties 
of the SSH model in the adiabatic limit by using the 
Kubo formula for the conductivity is quite important. 
We will show that, at the thermodynamical equilibrium 
for both electrons and lattice degrees of freedom and in- 
cluding all the vertex rinormalizations, the mobility is 
dominated, as expected jii by an "ad hoc" broadening 
of the energy levels. A proper choice of this broadening 
energy, taking into account in a qualitative way the miss- 
ing energy scales in the adiabatic limit, is able to recover 
the power-law observed in the experiments. Then, we 
discuss in some detail how the coupling strength among 
charge carriers and phonons is able to modify temper- 
ature dependence of the mobility and charge carrier lo- 
calization. We also measure the temperature dependent 
participation number showing a very weak temperature 
dependence in the experimental window and, on the con- 
trary, a very strong dependence on the coupling strength 
between charge carriers and intermolecular phonons. 

Furthermore, the analysis of the spectral function al- 
lows us to reconcile the results provided by the ARPES 
data (apparent band like description) with the computa- 
tional observation that the scattering length of the charge 
carrier becomes too short to be compatible with band 
transport £ 

Finally we focus on the optical conductivity of the 
model emphasizing the dependence on temperature and 
charge carrier density. This allows us to individuate a 
low energy peak, below any charge transfer excitation, 12 
that compares well with recent experimental results i 13 i 14 



II. THE MODEL 

The transport properties, the spectral function and 
optical conductivity of rubrene will be studied within 
the SSH model introduced by Troisi. 8 In this model the 
charge carriers move in a one-dimensional lattice hopping 
between next neighboring sites with a probability ampli- 
tude controlled by the relative position of the ions at the 
sites involved in the hopping. It can be summarized in 
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the following model hamiltonian: 

H = y E (^) 2 + \ E 0*) 2 + H el (1) {*})] = 

E(PA-^)KA|J|V)| 2 ^A-^+^) (4) 
where a, A' 

where 



= E M + a ( x i+l - x i)] ( c ! c *+l + c i+l c i) ■ ( 2 ) 



In eqs. (P) and (J2J) t is the bare electron hopping, 
a is the coupling constant that controls the link be- 
tween the electron hopping and ion displacement (xi) 
and, finally, m and k are the mass and the elastic con- 
stant of ions, respectively. We emphasize that the elec- 
tron dynamics is fully quantum (cj being the charge 
carrier creation operator) while the ion dynamics is as- 
sumed classic. The latter approximation is well justi- 
fied from the typical values of phonon frequencies ton 
and hopping constant t for rubrene. Following Ref.[16[ 
Hujq ~ 6meV and t ~ lAOmeV leading to an adiabatic ra- 
tio 7 = hu /t ~ 0.04. As we will discuss later, even if we 
are in a strong adiabatic regime the charge carrier mobil- 
ity can be still affected by very small quantum effects due 
to the one-dimensional nature of the model. On the other 
hand, the finite frequency conductivity and other proper- 
ties are not affected significantly by quantum fluctuations 
in this regime. In the following we will use dimensionless 
units measuring lengths in units of Iq = (fi/(2mcJo)) 1//2 
and energies in units of t. 

Taking advantage of the classic nature of the lattice 
distortions, the partition function can be written as 



Z = 



2mn 



L/2 



exp 



-/£Eo*) a 



Zel [{i 



(3) 

where Z e i[{xi}] is the quantum partition function of 
the electron subsystem given a deformation configuration 
{x^, L is lattice size and (3 = 1/KbT. 

Following Michielsen and de Raedtjil it is possible to 
estimate Z by using a Monte Carlo approach for the clas- 
sical degrees of freedom and exact diagonalization for 
the electron quantum dynamics. The method provides 
an approximation-free partition function of the model in 
the semiclassical limit. The only limitation is due to the 
computational time being controlled by the L x L ma- 
trix diagonalization. This constrains our analysis up to 
L = 128. In order to reduce the size effect, we use peri- 
odic boundary conditions. 

Within the same framework it is also possible to cal- 
culate the spectral function and the optical conductivity 
averaging the electronic properties at a given ion dis- 
placement configuration over the entire set of configura- 
tions weighted by Monte Carlo dynamics^ For instance, 
in the case of conductivity, for each configuration we cal- 
culate 



Px = 



exp[/3(£ A - m)] + 1 " 



(5) 



In eq.(dl) V = La is the system volume (a is the distance 
between next neighboring sites) and in eq.([S]) p, is the 
chemical potential, while E\ are the eigenvalues. 

The matrix element in eq. @ can be expressed in terms 
of the eigenstates of H e i. By using the unitary matrix 
that diagonalizes H e i, U(i,X), we can write 

(ai j\\>) = e u E [ u & + ^ *0 ( x \ c W' i y ) 

- tf(i + l,/i)tf(iy)(A|ct<v|A'}] (6) 
= E*i M*, A)I7(* + 1, A') - U(i + 1, X)U(i, A')] 



(7) 



where 

ti = -t + a (xi+i - Xi) . 
Then, the mobility can be defined as: 

Me = — linr Re\a(uj)] 
pe lj^o+ 



(8) 



where a{u) is obtained by averaging over the displace- 
ment configurations a(u; {xi}) and the density, p — 
N e /L, (N e being the charge carrier number) is 



(9) 



In the following we will focus our attention mainly to 
the limiting case in which a single electron is present in 
the system. ThenpA >-> exp[— j3E\]/Z e i in eq.(f4]) and the 
factor two, due to the spin degeneracy, drops out. 



III. TEMPERATURE DEPENDENCE OF 
MOBILITY 



As mentioned in the introduction, the mobility at the 
thermodynamic equilibrium, calculated assuming that 
the ion displacements are classical variables, is domi- 
nated by the energy broadening that we have to include 
in eq. ([?]). Indeed, in order to use that expression for any 
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finite lattice we have to replace the delta function with a 
Lorcntzian 



-d(E x - Ey +uj)^ — 77 

7T (E x - E x > + uy 



(10) 



The correct expression is obtained, then, for L i— > oo 
and 77 i — ^ . It can be easily shown that, as expected, 
this procedure leads to a vanishing mobility. Indeed, in 
the adiabatic limit, our calculation is equivalent to the 
classical problem of a particle in presence of off-diagonal 
disorder that is characterized by a vanishing mobility and 
localized eigen-functions except for that one correspond- 
ing to the zero energy eigenvalue! 19 ' 20 It is then clear 
that the system will exhibit a finite mobility only if we 
include a broadening representing the small (relevant) 
energy scale not included in the original model. In our 
opinion the most important missing energy scale is the 
quantum phonon ground state energy ftwo/2. The effect 
of this energy scale is usually negligible in the adiabatic 
limit, but becomes relevant in ID systems for the optical 
conductivity at energies less than Hujo/2 and, then, it is 
crucial for the mobility. 

In Fig.([T|) we show the temperature dependence of the 
mobility for the model parameters suitable for rubreneJ^ 
a/t = 0.09, u)o/t — 0.04. We restrict ourselves to the 
case of a single charge carrier and plot the mobility for 
different values of 77. It comes out that at 77 = Hujq/2 
a power-law behavior /i e ~ T~ s with 6 = 2.03 is re- 
covered. On the other hand, a smaller (larger) value of 
r\ provides a larger (smaller) exponent. In our opinion 
this is a strong indication that quantum fluctuations are 
indeed important in this system. It is also worth notic- 
ing that, at the thermodynamic equilibrium, the mobility 
does not depend on the ion mass [see eq.([3])]. Both the 
results, the need of an "external" energy scale and the 
mass independence of mobility, are not recovered within 
the dynamical approach proposed in Ref . Q pointing out 
that the dynamical and the equilibrium thermodynamic 
approaches are not equivalent. We also note that the 
mobility values we find are larger than those obtained in 
Ref.Q and surprisingly close to the experimental values. 

Summarizing, we have shown that the model of eq.([T]) 
in the adiabatic limit and at the thermodynamic equi- 
librium is able to recover the power-law observed in the 
experiments [see Fig. (JlJ] only introducing an "ad hoc" 
energy broadening that we associate to the quantum lat- 
tice fluctuations. On the other hand our analysis allows 
us to give a simple explanation of the physical mechanism 
responsible for the finite mobility. Due to the one dimen- 
sional nature of the model and the assumption that the 
lattice oscillations are classical variables all the charge 
carrier wave-functions are localized [see the snapshot re- 
ported in Fig. (J2J)]. In particular, as already proven in 
Ref.[l(|, a more detailed analysis shows that the wave 
functions whose energies are closer to band border are 
even more localized than those well within the energy 
band. The only exception to this description is given 
by the wave- function corresponding to vanishing energies 
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FIG. 1. Mobility vs Temperature for different values of r\. 
Mobility is in units of cm~ 2 /(s ■ Volt) and the temperature 
in Kelvin. The system size is L = 64. 
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n, 



<«i>r 

0.09 [ j 

0.06 

0.03 \ 
1 UmA 







a) 



0.6 r 

0.4 
0.2 




40 80 120 
C) 

E=-2.159 



40 80 120 




40 80 120 



FIG. 2. Snapshot of the electron wave-function at (typical) 
given deformation configuration at T = 347K. Panel a): av- 
erage density versus lattice sites; panel b)-d): average density 
at fixed electron energy versus lattice sites. The system size 
is L — 128 and the energy E is measured in units of t. 

This scenario leads to a vanishing mobility unless we 
assume an energy broadening that is able to provide the 
needed eigenvalue overlap. The power-law observed ex- 
perimentally in the temperature range 2fiw^ < KqT < 
5hujQ is, then, the result of entirely incoherent processes 
that have nothing to do with the band transport. It is 
worth noticing that this temperature behavior is due to 
the very small phonon energy involved in the scattering 
and it corresponds, in simple metals, to that typical of 
very high temperature regimei 21 ' 22 Finally, we empha- 
size that only a full quantum analysis will clarify how 
this scenario is modified at lower temperatures where the 
adiabatic approximation breaks down and the quantum 
fluctuations enter the problem in a more intrinsic way 
giving rise also to a Boltzmann-like contribution of the 
charge carriers that cannot be recovered in the present 
adiabatic approximation. 
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For a better understanding of the physical origin of the 
temperature dependence of the mobility we note that, as 
first observed by Troisi 8 and discussed in more detail by 
Fratini and Chichi^ the wave-function localization in- 
creases with temperature. In the equilibrium thermody- 
namic approach used in this paper the localization stems 
from the analogy with the disorder problem: increasing 
the temperature is equivalent to strengthen the disorder. 
On the other hand, if we measure the thermal average 
of any physical quantity the temperature enters not only 
through the distribution of ion displacements (as in the 
equivalent disorder problem), but also in the thermal av- 
erage of the physical quantity of interest at any fixed 
lattice configuration. Actually, with increasing the tem- 
perature, the thermal average involves more and more 
wave-functions corresponding to larger and larger eigen- 
values which correspond to less localized wave-functions. 
The competition between the two effects can give non 
trivial results. In order to clarify this issue we have cal- 
culated the temperature dependent participation number 
that provides a measure of the average effective localiza- 
tion: 



P = 



N 



(ii) 



P ranges from the number of the lattice sites (N) for 
delocalized states (traslational invariance) to 1 for a state 
fully localized on a single site. 

In Fig.© we report the temperature behavior of P 
showing that, in the temperature regime of interest, P 
does not exhibit a significant change. In this sense we 
think that the fact that the mobility decrease with the 
temperature is not strictly related to the wave-function 
localization, but rather it is due to an increase of the 
scattering rate among the electron and the lattice defor- 
mation. 
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FIG. 3. Participation number, P, as function of tempera- 
ture T for different values of the coupling constant a. The 
temperature, T, is measured in Kelvin and the system size is 
L = 64. 

On the other hand, it is interesting to emphasize the 
strong dependence of the participation number, P, with 
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the coupling constant a. In particular for a — 0.15 
P becomes equal to few lattice sites. At this value of 
the coupling constant the nature of the ground state has 
changed: the electron form a bond polaron32r— that is 
characterized by a very large effective mass. 

We end up this section presenting the mobility for dif- 
ferent coupling constants [Fig.©]. We observe that the 
power-law behavior is very robust and is recovered even 
in the most localized case reinforcing the idea that it is 
a purely incoherent mechanism. Finally we emphasize 
that, while the absolute value of the mobility decreases, 
the power 5 becomes smaller signaling a rather complex 
link between mobility and localization. 
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FIG. 4. mobility, [i e , as function of temperature T for differ- 
ent values of the coupling constant a. Mobility is in units of 
cm~ 2 /(s ■ Volt) and the temperature in Kelvin. The system 
size is L = 64. 



IV. SPECTRAL FUNCTION 

In the previous section we have shown that the equi- 
librium mobility in our model is due to purely incoherent 
processes and it cannot be ascribed to a simple band- 
like description. On the other hand, as mentioned in 
the introduction, ARPES measurements^ show that the 
effective electronic energy dispersion extracted from the 

dependent spectral function, A(k, oj), is very close to 
a simple cos(ka) band. This finding is usually considered 
as an indication in favor of a band-like scenario. In order 
to clarify this apparent contradiction and with the aim to 
validate the model studied in the present paper, we have 
calculated the spectral function for different k values and 
extracted the effective electronic band (see Fig.©). 

As it can be seen in the inset, the energy dispersion 
obtained following the main maxima of the spectral func- 
tion is very close to the bare band in agreement with the 
ARPES. However, a more careful analysis shows that the 
peaks exhibited by the spectral function are very broach 
making the quasi-particle description not well founded. 
Actually, for k — 7r/4 is evident a double peak struc- 
ture that persists for lower k values even if in a less ev- 
ident way. Finally, as expected, the spectral function 
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FIG. 5. The spectral function for different values of the 
wave-vector k. Measuring k in units of the inverse of the 
lattice spacing, a, we show from left to right the values 
ka = 0, ^, j. In the inset we report the effective elec- 
tron band (points) compared with the bare electronic band 
(full line). The temperature K B T/t = 0.21 and L = 64. 

presents an anomalous behavior at E = where the as- 
sociated wave-function is not localized] 19 ' 20 In our opin- 
ion, the present analysis provided a simple explanation of 
the apparent contradiction between ARPES and trans- 
port measurements: the absence of well defined quasi- 
particles makes the band-like description not applicable 
but, at the same time, the main peak of the spectral 
function inherits the bare band dispersion. Finally, we 
would mention that the effective energy dispersion de- 
pends very little with temperature in the range analyzed 
in this paper showing a very little increase of the band- 
width. In this sense the system is not characterized by 
any quasi-particle with heavy effective mass. 

V. OPTICAL CONDUCTIVITY 

As mentioned in the introduction, measurements of the 
optical conductivity (OC) are an important tool to inves- 
tigate the properties of OFET devices i 13 ' 14 For this rea- 
son we have calculated OC within the studied model. 
In Fig.© we show the OC for the parameter values 
relevant.— The OC exhibits a clear peak at low energies 
huj ~ 0.2t whose intensity decreases with the tempera- 
ture moving slightly towards high energies. It is worth 
noticing that, unlike the mobility, the peak position does 
not depend on the broadening energy that we still choose 
equal to Huiq/2. The result is of some interest since there 
are experimental evidenc o 13 ' 14 that, indeed, a peak is 
present at energies about 62meV (500cm _1 ) lower than 
any charge transfer process.— Assuming t ~ lAOmeV our 
estimate is a factor two lower than the measured value. 
Our estimation has to be considered quite reasonable for 
the very simple one-dimensional model we adopted. Ac- 
tually, as we will show in the following, our estimate can 
be even improved modifying slightly the phenomenolog- 
ical parameters of the model. 



We observe that the OC of the model exhibits also 
smaller structure at higher energies. However, for en- 
ergy very much larger than the bare charge carrier hop- 
ping, t, the model cannot be trusted in the framework of 
the rubrene OFET since at higher energies many charge 
transfer processes are observed 1 ^ that are not taken into 
account in the model presented here. Nevertheless these 
higher energies structures are still of interest for the 
model itself that has been proposed in many contexts. 




FIG. 6. The dimensional optical conductivity, E = 
a{u>)hV/[2n(ea) 2 ], as a function of the dimensionless pho- 
ton energy ui = Ml/t for different temperatures (measured 
in Kelvins). The lattice size is L — 64. 

As in the case of the mobility, it is interesting to study 
how the OC changes with the strength of the coupling 
constant, a, which, as discussed in the previous sec- 
tion, makes the system more and more "localized". As 
shown in Fig.©, when a increases, all the spectral weight 
move towards higher energies. In particular, at a — 0.12 
[Fig.©], the low energy peak move to w ~ 0.5t, becomes 
broader and looses intensity. This analysis suggests that 
a larger value of the low energy peak can be obtained by 
tuning the value of the charge-lattice coupling providing 
a better agreement with the experimental data. A fur- 
ther increase of a drives the system towards a localized 
state associated with a very large increase of the effec- 
tive mass and a clear optical gap opens up at low energies 
signaling the bond polaron formation (a = 0.15). 



VI. TRANSPORT PROPERTIES AND OPTICAL 
CONDUCTIVITY AT FINITE DENSITY 

Up to now we have focused our attention on the case 
of a single particle interacting with the lattice fluctua- 
tions, but the experiments, of course, are performed at 
a finite (even if small) particle density. Can we expect 
significant differences? The approach used in this work 
has the advantage that can deal very naturally with a 
finite number of particle and, then, we can address this 
point. In Fig.© we show the mobility as a function of the 
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FIG. 7. The dimensionless optical conductivity, E = 
a(u})hV/[2Tr(ea) 2 ], as a function of the dimensionless pho- 
ton energy ui — hQ,/t for different temperature (measured in 
Kelvins). The lattice size is L = 64. Upper panel a = 0.12, 
lower panel a — 0.15 



charge carrier density, p for low densities (p < 0.12) and 
different temperatures. We restrict our analysis to small 
density because OFET are characterized by extremely 
low density and also since we completely neglect charge- 
charge interaction in the model studied. It is clear that 
the mobility is significantly affected by the carrier density 
reducing its value when p increases. The single particle 
case studied in the previous sections provides, then, the 
maximum value for the mobility in this model. The non 
uniform behavior with temperature also signals that the 
temperature exponent of the mobility changes with the 
density. From our analysis stems out that, although the 
density increase provides an obvious increase of the con- 
ductivity (it is proportional to p), the mobility decrease 
compensates such increase and can even cancel it. There- 
fore, we expect a less linear increase of the conductivity 
with density. We interpret this result as the effect of 
the charge-charge effective interaction (mediated by the 
lattice fluctuations) that represents a further scattering 
mechanism for the charge carriers. The sensitivity of the 
mobility to the effective interaction suggests that the in- 
clusion of direct Coulomb interaction among the charges 
could be important even at low densities. 

We end up this section presenting the OC at differ- 
ent densities for the model parameters appropriate for 
rubrene. We still get a low energy peak as in the case 
of a single particle, but the peak position moves towards 
higher energies and the intensity decreases [Fig. ([5])] in- 
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FIG. 8. Mobility vs. charge carrier density, p for different 
temperature [T = 3487^ full line, T = 3UK dashed line, 
T = 2SIK dotted line, T = 2487f dashed-dotted line, T = 
215-ff long dashed line, T = 1&2K long-long dashed line]. The 
lattice size is L = 64. 



creasing the density p. This behavior, as it has been 
shown in Fig.©, is accompanied by a decrease of the 
mobility. We note that the observed behavior at finite 
density reminds that obtained increasing the coupling a 
[see Fig.©]. 
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FIG. 9. Optical conductivity, £ = a(uj)hV/[2n(ea) 2 ], vs. di- 
mensionless photon energy, uj — hQ/t for different charge car- 
rier density p [p = 0.014 full line, p — 0.014 dashed line, 
p = 0.014 dotted line, p = 0.014 dashed-dotted line]. The 
temperature is T = 3487^. Energy are measured in units of 
the bare hopping, t. 



VII. CONCLUSIONS 

In this work we have studied the thermodynamic equi- 
librium properties of the SSH model in ID assuming 
the ion displacements as classical variables. We focused 
our attention on mobility, spectral function and optical 
conductivity in the adiabatic regime (7 = 0.04). Actu- 
ally, this regime is particularly interesting for its con- 
nection with the low energy physics of OFET based on 
Rubrene single crystal as proposed in Refs|H andHiil. We 
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find that, as expected at the thermodynamic equilibrium, 
the model mobility vanishes unless an "ad hoc" energy 
broadening is introduced in the model. Interestingly if we 
choose this energy broadening of the order of Hujq/2 we 
are able to recover the temperature dependence observed 
in the experiments. This result suggests that the lattice 
quantum fluctuations not considered in our approach are 
crucial for the mobility even in a strong adiabatic regime. 
It is worth noticing that the present approach shows a 
significant difference with the approximated dynamical 
approach£ii£ where a finite mobility is obtained without 
invoking an "ad hoc" energy broadening. In our opin- 
ion the difference stems from the fact that the dynamics 
adopted does not bring the system to the full thermody- 
namic equilibrium. 

We also showed that the single-particle spectral func- 
tion, A(k,ui) exhibits very broad peaks that cannot be 
described as a simple quasi-particle. However, the posi- 
tion of the main peak inherits the bare electron energy 
dispersion. This result could explain the apparent con- 
tradiction between ARPES and the power law tempera- 
ture of the mobility observed in the transport measure- 
ment of Rubrene OFET. 



Finally, the analysis of the OC of the model suggests 
the existence of a low energy peak in the energy range 
0.2t < ui < 0.6t depending on the charge-lattice cou- 
pling a and the charge carrier density. As discussed in 
the previous sections, this result is of some interest since 
there are experimental evidence on single crystal Rubrene 
OFE T 13 i 14 showing that, indeed, a peak is present at 
energies about 62meV lower than any charge transfer 
process Assuming in our model a bare charge hopping, 
t ~ 140meV, our estimate is in a reasonable agreement 
with the measured value. More experimental data, in 
particular temperature dependence of the peak, and more 
realistic models will be needed for a full characterization 
of this low energy OC peak. 
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